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摘要 : 为 提高 对 流 层 天 顶 延 迟 (Zenith Tropospheric Delay，ZTD ) 估计 精度 ， 基 于 传 
统 对 流 层 天 顶 延 迟 建 模 思 路 ， 提 出 Saastamoinen, Askne 和 GPT3 多 模型 融合 的 对 流 层 天 顶 
延迟 估计 方法 。 分 别 采 用 Saastamoinen 和 Askne 模型 估计 干 延 迟 和 湿 延 迟 ， 并 引入 GPT3 
模型 提供 温度 、 气 压 、 水 汽 压 、 大 气 加 权 平 均 温度 和 水 汽 垂 直 递 减 率 气 象 参数 。 利 用 全 球 大 
地 测量 观测 系统 (Global Geodetic Observing System, GGOS ) Atmosphere 和 国际 GNSS 服 
务 机 构 (International GNSS Service，IGS ) 提供 的 亚洲 区 域 66 个 IGS 站 的 2016-2018 年 对 
流 层 天 顶 延 迟 数据 对 本 文 方法 进行 评估 ， 结 果 表明 ， 以 GGOS Atmospheres 数据 为 参考 时 ， 
Sas+Ask+GPT3 模型 精度 (RMS: 4.53cm ) 较 同 等 条 件 下 的 Sas+Ask+UNB3m 和 Sas+GPT3 
模型 分 别提 高 约 29% 和 19%, VAIGS 对 流 层 天 顶 延迟 数据 为 参考 时 ，Sas+Ask+GPT3 模型 
精度 (RMS: 4.35cm ) 较 另 两 种 模型 分 别提 高 约 25% 和 14%; SastAsk+GPT3 模型 误差 具有 
冬季 小 而 夏季 大 的 季节 性 特征 ， 且 夏季 误差 明显 低 于 另 两 种 模型 ; 在 空间 分 布 上 ， 模 型 误差 
随 高 程 或 纬度 的 增加 而 降低 。 研究 表明 ,本 文 方法 可 用 于 亚洲 区 域 对 流 层 天 顶 延 迟 的 高 精度 
实时 修正 。 
关键 词 : 对 流 层 天 顶 延 迟 ; GPT3 模型 ; UNB3m 模型 ; Saastamoinen 模型 ; Askne 模型 


中 图 分 类 号 : P228.4 文献 标识 码 : A 文章 编号 : 1672-7673 (2023) 01 


由 于 地 球 周围 的 对 流 层 对 电磁 波 的 折射 效应 ,电磁波 信号 的 传播 速度 变 慢 ,传播 路 径 齐 
曲 , 导致 电磁 波 信号 传播 延迟 , 这 种 现象 称 为 对 流 层 延 迟 .对流层 延迟 可 分 为 静 力 学 延迟 (也 
称 对 流 层 干 延 迟 ，Zenith Hydrostatic Delay (ZHD ) ) 与 非 静 力 学 延迟 (又 称 对 流 层 湿 延 迟 ， 
Zenith Wet Delay (ZWD) ) ， 其 中 ， 对 流 层 干 延迟 由 大 气 中 的 干燥 气体 引起 ， 占 总 延迟 的 
90 多 左右， 变化 比较 稳定 ; 对 流 层 湿 延迟 由 水 汽 引 起 ， 占 总 延迟 的 10% 左 右 ， 具 有 时 变性 。 


在 对 流 层 延迟 改正 中 ， 通 常 将 电磁 波 信号 在 天 顶 方 向 的 延迟 通过 映射 函数 投影 至 传播 路 径 ， 
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为 2m， 在 地 平方 向 可 达 20ml1， 因 此 ， 对 流 层 延迟 是 全 球 卫星 导航 系统 (Global Navigation 
Satellite System, GNSS) 和 其 长 基线 干涉 测量 (Very Long Baseline Interferometry, VLBI) 
等 空间 大 地 测量 技术 的 重要 误差 源 之 一 上 1。 

对 流 层 是 非 色散 介质 ， 无 法 采用 频率 组 合 的 方式 消除 对 流 层 延迟 ,通常 用 对 流 层 延迟 经 
验 模型 予以 修正 。 当 前 常用 的 对 流 层 天 顶 延迟 经 验 模 型 大 致 可 以 分 为 两 类 , 一 类 是 基于 实测 
气象 参数 的 Saastamoinen 和 Hopfield 等 模型 ， 另 一 类 是 仅 需 测 站 地 理 位 置信 息 的 非 气 象 参 
数 模 型 ， 如 美国 广 域 增强 系统 (Wide Area Augmentation System, WAAS ) 采用 的 UNB 模型 
多 和 欧洲 广 域 增强 系统 采用 的 EGNOS RANI, RFK MARSA MM AAA THER RA! 44 
改正 效果 高 度 依赖 于 实测 地 表 气 象 参 数 ， 利 用 实测 地 表 气 象 参 数 估计 对 流 层 天 项 延迟 ,改正 
精度 可 达 厘 米 级 外， 车 利用 标准 气象 参数 则 对 流 层 天 顶 延 迟 改正 效果 较 差 ， 这 限制 了 此 类 模 
型 在 实时 定位 、 和 导航 与 授时 中 的 推广 应 用 ， 并 且 Saastamoinen 和 Hopfield 等 模型 对 对 流 层 
天 顶 延 迟 处 理 比较 简单 ， 使 得 对 流 层 天 顶 延 迟 改正 精度 难以 有 较 大 提高 ; UNB 系列 模型 和 
EGNOS 模型 无 需 实 测 地 表 气 数 参 数 ， 使 用 简便 ， 它 们 的 全 球 平均 对 流 层 天 顶 延 迟 改正 精度 
可 以 与 基于 实测 气象 参数 的 Saastamoinen 和 Hopfield 模型 相当 ， 然 而 ， 由 于 这 些 模型 将 大 
地 按 15" 的 纬度 间隔 划分 ， 分 状 率 较 差 ， 且 忽略 经 度 的 影响 ， 导 致 这 些 模型 无 法 细致 地 刻画 
对 流 层 延迟 的 局 部 变化 ， 在 某 些 区 域 的 改正 效果 不 佳 "。 

考虑 到 现 有 对 流 层 天 顶 延 迟 经 验 模 型 的 局 限 性 , 采用 与 实测 气象 参数 相符 合 的 气象 参数 
模型 估计 对 流 层 天 顶 延迟 成 为 改善 对 流 层 天 顶 延迟 改正 效果 的 一 种 有 效 方法 , 在 近年 来 得 到 
1 ERR. x [8148 GPT 与 UNB3m 气象 参数 模型 估计 的 气温 、 大 气压 和 水 汽 压 参 
数 作 为 Saastamoinen 模型 的 输入 , 有 效 提 高 了 Saastamoinen 模型 的 适用 性 ; 文 [1] 利用 GPT2 
模型 提供 的 气温 、 大 气压 和 相对 湿度 应 用 于 Saastamoinen 模型 中 以 估计 对 流 层 天 顶 延 认 ， 
建立 了 适用 于 全 球 范 围 的 对 流 层 天 顶 延 迟 精 化 模型 ; x [9-10] 通 过 GPT2w 模型 计算 气象 参 
数 ， 并 分 别 将 气象 参数 作为 Hopfield 和 Saastamoinen 模型 的 输入 ， 摆 脱 了 Saastamoinen 和 
Hopfield 模型 对 实测 气象 参数 的 依赖 , 提高 了 两 种 经 典 模型 的 普 适 性 。 总 结 已 有 的 研究 发 现 ， 
通过 GPT/GPT2/GPT2 和 UNB3m 模型 能 够 提高 Hopfield 和 Saastamoinen 模型 的 适用 性 , 但 
由 于 它们 的 湿 延 迟 估 计 精 度 较 差 , 这 限制 了 两 种 经 典 模型 对 流 层 天 顶 延 迟 估 计 精 度 的 进一步 
提高 ; 此 外 ， 不 同 气 象 参数 模型 在 不 同 区 域 都 具有 一 定 适 用 性 。 

为 进一步 提高 基于 气象 参数 的 对 流 层 天 顶 延迟 模型 的 改正 精度 , 本 文 在 上 述 研 究 的 基础 


上 ， 将 目前 最 新 的 GPT 系列 模型 一 GPT3 模型 提供 的 气象 参数 应 用 于 对 流 层 天 顶 延迟 计算 ， 


其 中 ， 采 用 Saastamoinen 模型 估计 对 流 层 干 延 迟 ， 采 用 性 能 更 好 的 Askne 模型 估计 对 流 层 
湿 延 迟 ，Saastamoinen 模型 所 需 的 气压 和 Askne 模型 所 需 的 水 汽 压 、 大 气 加 权 平 均 温度 与 
水 汽 垂直 递减 率 4 个 关键 气象 参数 均 由 1 x 1° 空间 分 辨 率 的 GPT3 模型 提供 ， 由 此 构成 
Sas+Ask+GPT3 融合 模型 .考虑 到 GPT3 模型 在 亚洲 区 域 的 对 流 层 天 顶 延迟 估计 中 应 用 较 少 ， 
以 GGOS Atmosphere fe IGS 发 布 的 对 流 层 天 顶 延 迟 产 品 作 为 参考 ， 分 析 评 估 了 组 合 模 型 和 
同等 条 件 下 的 SastGPT3 和 SastAsk+UNB3m 应 用 于 亚洲 区 域 的 对 流 层 天 顶 延迟 改正 精度 ， 
并 分 析 了 模型 误差 的 时 空 分 布 特征 。 研 究 结 果 表 明 ，GPT3 模型 提供 的 气象 参数 相对 更 符合 

亚洲 区 域 的 实际 情况 ，Sas+Ask+GPT3 融合 模型 对 流 层 天 顶 延 迟 的 估计 精度 总 体 最 优 ， 可 应 
用 于 亚洲 区 域 的 对 流 层 天 顶 延迟 的 实时 改正 ， 并 有 具有 所 需 参 数 少 和 使 用 方便 的 优点 。 
1 多 模型 融合 对 流 层 天 顶 延 迟 估计 方法 

Saastamoinen 和 Hopfield 模型 是 两 种 应 用 广泛 的 对 流 层 天 顶 延 迟 经 验 模型 ， 其中， 前 者 
精度 受 高 程 影响 很 小 ， 而 后 者 精度 随 高 程 的 增加 而 降低 ， 故 Saastamoinen 模型 的 适用 范围 
3-9, [e Saastamoinen 模型 对 对 流 层 湿 延 迟 处 理 较为 简单 ， 估 计 误 差 较 大 。Askne 模型 融 
入 大 气 加 权 平 均 温度 和 水 汽 垂直 递减 率 参 数 0 1， 对 流 层 湿 延 迟 估计 精度 更 高 ， 但 这 些 参数 
依赖 于 探 空气 球 测量 , 限制 了 其 在 对 流 层 天 顶 延 迟 实时 改正 中 的 应 用 。 GPT3 模型 在 GPT2w 
模型 基础 上 又 引入 了 梯度 算法 ， 进 一 步 优 化 了 模型 性 能 ， 可 以 提供 全 球 范围 内 1° x 1° 高 分 
辨 率 的 气象 参数 值 04。 为 提高 对 流 层 天 顶 延 迟 经 验 模 型 的 估计 精度 和 适应 能 力 ， 本 文 提出 
一 种 多 模型 融合 的 对 流 层 天 顶 延 迟 估 计 方 法 ， 分 别 采 用 Saastamoinen 和 Askne 模型 估计 对 
流 层 干 延迟 和 对 流 层 湿 延 迟 ， 两 种 模型 所 需 的 气象 参数 全 部 由 GPT3 模型 提供 ， 在 提高 对 流 
层 天 顶 延 迟 估计 精度 的 同时 实现 对 流 层 天 顶 延 迟 的 实时 改正 。 
1.1 Saastamoinen 模型 

Saastamoinen 模型 将 对 流 层 分 为 两 层 进行 积分 , 第 1 层 为 从 地 表 到 10 km 高 度 的 对 流 层 
顶 ， 该 层 的 温度 变化 率 为 6.5'C /km; 第 2 层 为 从 10 km 到 50 km 高 度 的 平流 层 顶 ， 将 该 层 
的 温度 视 为 常数 。Saastamoinen 模型 首次 将 被 积 子 数 按照 天 顶 距 三 角 函 数 展 开 逐 项 进行 积 
分 ， 并 把 对 流 层 天 顶 延 迟 分 为 对 流 层 干 延 迟 和 对 流 层 湿 延 迟 两 项 分 量 之 和 ,两 项 分 量 的 表达 
A) 


ZHD = DIEM 


floh) 
ZWD = 0.002277. (0.05 + 128) ; (1) 
fly, h) = An :0.0020600(2/) 0.00028 


KB, P, T, €, PARAR ZB AJ& (hPa) . WAI (K). KA (hPa) 、 测 站 


纬度 (rad) 和 高 程 (km ) 。 
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1.20 Askne 模型 
Askne 模型 融入 大 气 加 权 平 均 温 度 和 水 汽 重 直 递减 率 参 数 以 改善 对 流 层 湿 延 迟 估计 精 
度 ， 对 流 层 湿 延 迟 表 达 式 为 上 


=). — (2) 


ZWD = 107 (k; + ; 


(A+ 1)9m 
其 中 ，Tim 为 大 气 加 权 平 均 温 度 (K); AJ KAE AMRF (hPa/m ) ; 9m 为 重力 加 速度 ， 
gm = 9.80665 m/s?; ks = 16.52K/hPa; ks = 3.776 x 10°K?/hPa; Rg = 287.0538 J/(K - g), 
1.3 GPT3 模型 
GPT 系列 模型 是 Boehm 等 利用 欧洲 中 尺度 天 气 预 报 中 心 (European Centre for 
Medium-Range Weather Forecasts, ECMWF) 长 期 的 再 分 析 气 象 资料 建立 的 全 球 气 象 参数 经 
验 模型 , 仅 需 知道 测 站 地 理 位 置信 息 与 年 积 日 便 可 以 获得 地 表 温度 、 大 气压 力 和 水 汽 压 等 气 
象 参 数 , 在 全 球 范围 内 得 到 广泛 应 用 。 GPT3 模型 是 GPT 系列 模型 的 最 新 版 本 , 它 在 GPT2w 
模型 的 基础 上 引入 了 梯度 算法 以 进一步 优化 模型 ， 可 以 提供 1° x 1° 空间 分 辨 率 的 格 网 气象 
参数 值 。 格 网 点 气象 参数 的 计算 公式 为 
doy 


doy doy doy 
27) + B,sit 27) + Ascos( - 47) + Bosin( 4m) 
365.257) * Pisln( 365 25 cos (565.25 7) + Basin(365.2547) (3) 


其 中 ，E 为 格 网 点 气象 参数 ; doy 为 年 积 日 ; A 为 气象 参数 的 平均 值 ，A1 和 BI 为 气象 参数 的 
年 周期 振幅 ;42 和 2 为 气象 参数 的 半年 周期 振幅 。 
1.4 融合 模型 

针对 Saastamoinen 和 Askne 模型 对 实测 气象 参数 的 依赖 性 ， 本 文 提出 一 种 适用 于 亚洲 
区 域 的 对 流 层 天 顶 延迟 融合 估计 模型 ,该 模型 融合 GPT3 模型 提供 的 大 气压 力 、 地 表 温 度 和 
水 汽 压 等 高 精度 气象 参数 来 优化 对 流 层 干 延 迟 和 对 流 层 湿 延 迟 估计 中 所 需 的 气象 参数 .从 而 
克服 对 实测 气象 参数 的 依赖 ,提高 对 流 层 干 延 迟 和 对 流 层 湿 延 迟 估计 精度 ， 同 时 实现 对 流 层 


天 顶 延 迟 的 实时 改正 。 融 合 模型 的 对 流 层 天 顶 延 迟 计算 步 骤 为 (1) 根据 测 站 经 度 、 纬 度 


= Ay + 4 A; cos( 


高 程 和 年 积 日 信息 计算 距离 测 站 邻近 的 4 个 格 网 点 的 气压 、 温 度 、 水 汽 压 、 水 汽 垂直 递减 率 
和 大 气 加 权 平 均 温 度 气 象 参数 ， 并 利用 (4) 一 (7) 式 将 格 网 点 上 的 气压 、 温 度 、 水 汽 压 进 
高 程 修正 ， 将 其 归 算 为 测 站 高 程 上 的 格 网 点 参数 ; (2) 将 测 站 高 程 上 的 4 个 格 网 点 参数 
的 气象 参数 进行 双 线 性 插值 ， 获 得 测 站 位 置 的 气象 参数 ; (3) 将 插值 得 到 的 气象 参数 带 入 
模型 来 计算 对 流 层 干 延 迟 和 对 流 层 湿 延 迟 。 


P = Py-exp(—gm -dMir/(R, - T,))/100 
T, = Ty - (1 + 0.6077Q) 
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T = To 4 dT -dh (5) 
e = Q- P/(0.622 + 0.378Q) (6) 
其 中 ， 硬 和 T0 分 别 为 格 网 点 上 的 温度 (K) 和 气压 (hPa) ; dT 为 温度 垂直 递减 率 ; CA 
X; 有 为 高 程 修正 量 (m); We 和 证 分 别 为 大 气 摩 尔 质 量 与 气体 常数 ， 
dM, = 0.028965 kg/mol, Ry =8.3143J/(K mol， 其 他 参数 含义 同上 文 。 
2 精度 评定 方法 
2.1 研究 区 域 
亚洲 地 域 东 西南 北 跨度 大 ,气候 复杂 多 变 ， 为 验证 融合 模型 在 亚洲 区 域 的 对 流 层 天 顶 延 
迟 改 正 有 效 性 ， 本 文选 取 亚 洲 区 域 66 IGS 站 对 对 流 层 天 顶 延迟 改正 效果 进行 分 析 评 估 研 
jb. 这 些 测 站 具备 地 域 跨度 大 、 海 拔 与 气候 差异 明显 等 特点 ,研究 区 域 的 范围 为 10°N-63°N， 
66"E-159*E， 选 取 的 亚洲 区 域 66 个 IGS 站 的 地 理 分 布 如 图 1, SIGS 站 的 经 纬度 和 高 程 数 


值 可 通过 网 址 https://igs.org/network/#new-site-checklist 查询 ， 限 于 篇 幅 ， 在 此 不 再 列 出 。 


Al 亚洲 区 域 1GS 站 点 分 布 


Fig.l Distribution of the selected IGS stations over Asia 
2.2 数据 来 源 
GGOS Atmosphere 和 IGS 提供 的 对 流 层 天 顶 延 迟 产品 具有 很 高 的 质量 ， 可 以 作为 对 流 
层 天 顶 延 迟 研 究 与 应 用 的 标准 数据 .本 文 利用 这 两 种 对 流 层 天 顶 延 迟 产 品 作为 参考 值 对 融合 
模型 的 精度 进行 分 析 评估 ， 其 中 ，GGOS Atmosphere 基于 ECMWF 40a 再 分 析 资料 ， 可 以 提 
供 1979 年 至 今 空 间 分 辨 率 为 2.5°x2° ( 经 度 x 纬 度 ) 、 时 间 分 辨 率 为 6h (每 天 0: 00、6: 00, 


12: 00 和 18: 00) 的 全 球 对 流 层 天 顶 延 迟 格 网 数据 ， 数 据 可 以 从 https:Wvmf.geo.tuwien.ac.at 


获取 ; IGS 提供 的 全 球 对 流 层 天 顶 延 迟 数据 由 GNSS 解 算 ， 时 间 分 辩 率 为 Smin， 数 据 可 以 
从 https://cddis.nasa.gov 获取 。 本 文选 取 2016-2018 年 的 对 流 层 天 顶 延 迟 产 品 作 为 参考 , 分 析 
评估 了 融合 对 流 层 天 顶 延 迟 模型 误差 的 时 空 分 布 特征 。 
2.2 精度 评定 指标 

分 别 以 GGOS Atmosphere Fe IGS RAP ART AAA HAIER BABA, ARP 
& (bias) 与 均 方 根 误差 (RMS) 两 个 指标 对 不 同 模型 的 对 流 层 天 顶 延 迟 估计 结果 进行 分 析 


评估 ， 其 计算 公式 分 别 为 


4 


bias = 一 (ZTD? — ZTD?) ; (7) 


N 
1 sesamin. An 
RMS |e ZTD?)? ; (8) 


其 中 ，ZTD? 为 模型 计算 值 ; ZTD? 为 GGOS Atmosphere 和 IGS 提供 的 对 流 层 天 顶 延 迟 数 据 ; 
六 为 对 流 层 天 项 延迟 数据 的 样本 个 数 。 
3 模型 验证 

为 更 好 的 分 析 评 估 Sas+Ask+GPT3 融合 模型 在 亚洲 区 域 的 对 流 层 天 项 延迟 改正 的 有 效 
性 ， 在 国定 GPT3 模型 的 条 件 下 与 Saastamoinen 模型 形成 Sas-GPT3 模型 ， 其 中 。 干 、 湿 延 
iR 34 4]Jf] Saastamoinen 模型 计算 ， 所 需 参 数 由 GPT3 模型 提供 ; 在 固定 SastAsk 模型 的 条 件 
下 与 UNB3m 模型 形成 SastAsk+UNB3m 模型 ， 其 中 ， 干 、 湿 延迟 分 别 利用 Saastamoinen 
和 Askne 模型 计算 ， 所 需 参 数 由 UNB3m 模型 提供 ， 并 将 本 文 融合 模型 与 另外 两 种 模型 的 对 
流 层 天 顶 延迟 估计 误差 进行 对 比分 析 。 
3.1 平均 偏差 和 均 方 根 误差 统计 结果 

分 别 以 GGOS Atmosphere 和 IGS 发 布 的 2016-2018 年 对 流 层 天 顶 延 迟 数据 作为 参考 ， 
计算 3 种 融合 模型 在 66 个 IGS 站 2016-2018 年 的 平均 偏差 与 均 方 根 误差 ， 如 图 2 和 图 3， 
66 个 IGS 站 的 平均 偏差 与 均 方 根 误差 统计 结果 见 表 1。 

由 表 1 可 知 ， 以 GGOS Atmosphere 产品 为 参考 时 ，Sas+Ask+GPT3 模型 的 平均 偏差 相 
对 同等 条 件 下 的 SastAsk+UNB3m 和 Sas+GPT3 模型 分 别 降 低 约 41% 和 72%, 均 方 根 误 差 相 
对 另 两 种 模型 分 别 降 低 约 29% 和 19%; 以 IGS 对 流 层 天 顶 延 迟 产品 为 参考 时 ， 
Sas+Ask+GPT3 模型 的 均 方 根 误差 相对 另 两 种 模型 分 别 降低 约 25% 和 14%， 平 均 偏差 相对 
Sas+GPT3 模型 降低 约 78%， 但 其 平均 偏差 大 于 同等 条 件 下 的 Sas+AsktUNB3m 模型 。 由 图 


2(a) 和 图 3(a) 可 以 看 到 ， 无 论 是 以 GGOS Atmosphere 还 是 IGS 对 流 层 天 顶 延 迟 数据 为 标准 ， 


3 种 模型 的 平均 偏差 均 存在 很 多 负 值 ， 导 致 66 个 IGS 站 的 平均 偏差 也 为 负 值 ; 
Sas+Ask+UNB3m 模型 的 平均 偏差 为 负 值 的 情况 少 于 Sas+Ask+GPT3 和 Sas+GPT3 模型 ， 特 
别 是 在 以 IGS 对 流 层 天 顶 延 迟 产 品 为 参考 时 ， 该 模型 在 39 个 IGS 站 的 平均 偏差 为 负 值 , 在 
27 个 站 为 正 值 ， 而 另 两 种 模型 在 大 部 分 IGS 站 为 负 值 ， 导 致 根据 式 (7) 计算 的 
Sas+Ask+UNB3m 模型 的 平均 偏差 偏 小 ; 图 2 和 图 3 还 显示 3 种 模型 在 局 部 水 汽 丰富 且 变 化 
复杂 的 区 均 会 出 现 较 大 偏差 (如 HKSL.WUHN 和 SHAO 等 IGS 附近 区 域 ), 但 Sas+Ask+GPT3 
模型 的 精度 相对 最 优 ， 说 明 Sas+Ask+GPT3 模型 效果 明显 。 

Sas+Ask+GPT3 和 Sas+Ask+UNB3m 模型 相当 于 在 国定 对 流 层 天 顶 延 迟 模型 的 条 件 下 与 
不 同 气象 参数 模型 形成 的 融合 模型 ,而 Sas+Ask+GPT3 和 Sas+GPT3 模型 相当 于 在 固定 气 
参数 模型 的 条 件 下 与 不 同 对 流 层 天 顶 延 迟 模型 形成 的 融合 模型 , 故 Sas+Ask+GPT3 模型 是 在 
同等 条 件 下 与 另 两 种 模型 进行 对 比 。Sas+Ask+GPT3 模型 的 对 流 层 天 顶 延迟 改正 效果 显著 优 
于 另 两 种 模型 ， 这 一 方面 说 明 Askne-ZWD 模型 的 准确 性 比 Saastamoinen-ZWD 模型 高 ， 另 
一 方面 说 明 GPT3 模型 在 亚洲 地 区 比 UNB3m 模型 具有 更 好 的 适应 性 ， 这 是 因为 UNB3m 模 
型 利用 北美 地 区 的 探 空 数据 构建 而 具有 一 定 的 区 域 性 。 此 外 ， 整 体 而 言 ，Sas+GPT3 模型 精 
度 高 于 Sas+Ask+UNB3m 模型 ,说明 对 流 层 天 顶 延 迟 估计 精度 不 仅 和 对 流 层 天 顶 延 迟 模型 的 
准确 性 有 关 ， 还 取决 于 气象 参数 的 准确 性 ， 间 接 验 证 了 GPT3 模型 计算 气象 参数 的 有 效 性 。 
本 文 将 对 3 种 模型 对 流 层 天 顶 延 迟 估计 误差 的 时 空 变化 特征 进行 分 析 。 

表 1 3 种 融合 模型 的 对 流 层 天 顶 延 迟 估计 误差 统计 值 (单位 : cm) 


Table 1 Error statistics of the ZTD estimates derived three combined models (unit: cm) 


Sas+Ask+GPT3 Sas+Ask+UNB3m Sas+GPT3 
Reference data 
bias RMS bias RMS bias RMS 
GGOS Atmosphere -0.61 4.35 -1.04 6.09 -2.16 5.37 


IGS -0.35 4.53 -0.20 6.03 -1.59 5.27 
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图 2 相对 GGOS Atmosphere ;* 49%} LACK TRAE IR ed FI hE 5 3] Z RR E 


Fig.2 Bias and RMS between the ZTD estimates and GGOS Atmosphere products 


图 3 相对 IGS 产品 的 对 流 层 天 顶 延 迟 估计 平均 偏差 与 均 方 根 误差 


Fig.3 Bias and RMS between the ZTD estimates and IGS products 


32 平均 偏差 与 均 方 根 误差 的 时 间 分 布 特征 
ADM 3 种 模型 平均 偏差 与 均 方 根 误差 每 日 的 变化 情况 ， 以 IGS 对 流 层 天 顶 延 迟 数据 
为 参考 ， 将 2016-2018 年 66 个 IGS 站 的 平均 偏差 和 均 方 根 误差 月 均值 进行 统计 ， 如 图 4。 
从 图 4 可 以 看 到 ，Sas+Ask+GPT3、Sas+Ask+UNB3m 和 Sas+GPT3 模型 的 月 均 平均 偏 
差 和 均 方 根 误 差 均 表现 出 显著 的 季节 性 变化 ， 月 均 平均 偏差 和 均 方 根 误差 在 夏季 7 月 或 8 
月 达到 最 大 值 ， 在 冬季 12 月 或 1 月 达到 最 小 值 ， 且 月 均 平均 偏差 和 均 方 根 误 差 在 夏季 变化 
度 相 对 较 大 ， 冬 季 变 化 较 小 ， 这 是 由 于 夏季 空气 湿度 较 大 ， 对 流 层 湿 延 迟 变 化 剧烈 ， 冬 季 


空气 比较 干燥 , 对 流 层 湿 延 迟 变化 较 小 ,而 对 流 层 湿 延迟 又 是 对 流 层 天 顶 延 迟 估计 误差 的 主 
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项 ; SastAsk+UNB3m 模型 的 月 均 平均 偏差 和 均 方 根 误差 在 3 种 模型 中 最 大 ， 特 别 是 在 水 汽 
变化 剧烈 的 夏季 ,表明 UNB3m 模型 预 估 的 气象 参数 出 现 较 大 偏差 ; 与 另 两 种 模型 相 比 较 而 
言 , SastAsk+GPT3 模型 的 月 均 平均 偏差 和 均 方 根 误差 相对 较 小 , 虽然 在 夏季 月 份 误差 增 大 ， 
但 各 月 份 之 间 变 化 相对 稳定 ， 说 明 GPT3 模型 预 估 的 气象 参数 比 UNB3m 模型 更 加 准确 ， 同 


时 说 明 即 使 在 水 汽 变 化 复杂 的 夏季 Askne 模型 也 能 较 好 地 反映 对 流 层 湿 延迟 变化 情况 。 


图 4 亚洲 区 域 IGS 站 ZTD 估计 月 均 平均 偏差 与 均 方 根 误差 

Fig.4 Monthly averaged bias and RMS of the ZTD estimates in the selected IGS stations over Asian 
3.3 平均 偏差 与 均 方 根 误差 的 空间 分 布 特征 
3.3.1 平均 偏差 与 均 方 根 误差 的 高 程 分 布 特征 

亚洲 区 域 地 势 高 、 地 表 起 伏 大 ， 对流 层 天 顶 延迟 变化 和 高 程 密切 相关 。 为 分 析 3 种 模型 
平均 偏差 与 均 方 根 误 差 在 高 程 上 的 变化 情况 , 将 66 35,505 3 £24 0 ~ 500m. 500 ~ 1000m, 
1500 ~ 2000m 和 >2000m 进行 划分 ， 对 2016-2018 年 各 范围 站 点 的 平均 平均 偏差 与 均 方 根 误 
差 进行 统计 ， 如 图 5， 由 于 在 亚洲 区 域 1000 ~ 2500m 高 程 范围 内 无 IGS 站 ， 故 未 统计 该 高 
程 范围 内 的 误差 。 

由 图 5 可 知 ，3 种 模型 的 平均 偏差 与 高 程 并 无 明显 关联 ， 但 在 高 程 大 于 2000m 时 
Sas+Ask+UNB3m 模型 的 平均 偏差 较 大 ， 可 能 是 在 该 高 程 范围 内 统计 的 站 点 个 数 很 少 的 原 
因 ; 相反 ，Sas+Ask+GPT3 和 Sas+GPT3 模型 的 平均 均 方 根 误差 与 高 程 有 着 明显 规律 ， 即 随 
着 高 程 的 增 大 ， 两 种 模型 的 均 方 根 误差 呈现 明显 的 递减 趋势 , 这 是 因为 高 海拔 地 区 的 折射 率 
的 积分 区 间 较 短 , 对流层 天 顶 延 迟 相 对 较 小 ; Sas+Ask+UNB3m 模型 在 高 程 ~ 500m 范围 内 
的 均 方 根 误差 较 大 (超过 6cm ) ， 当 高 程 大 于 500m 时 均 方 根 误差 明显 变 小 (小 于 Scm ) ， 
有 递减 趋势 ， 但 在 高 程 大 于 2000m 时 均 方 根 误差 又 突然 增 大 ， 同 样 可 能 是 因为 此 范围 内 站 
点 很 少 。 另 外 ， 从 图 5(a) 可 以 看 出 ， 当 高 程 小 于 500m 时 3 种 模型 的 平均 偏差 均 为 负 值 ， 说 


ERARA DURS ST MLAB ATER. ABP je, SastAsk+GPT3 模型 在 各 高 程 上 
估计 精度 均 优 于 另 两 种 模型 ， 特 别 是 对 于 2000m 以 上 的 高 程 其 均 方 根 误差 优 于 2cm, XX 
优 于 同等 条 件 下 的 Sast+Ask+UNB3m 模型 。 对 流 层 天 顶 延 迟 模型 易 受 高 程 影响 的 原因 在 于 ， 
对 流 层 天 顶 延 迟 分 布 和 高 度 密切 相关 ， 且 GPT3 和 UNB3m 气象 参数 模型 的 高 程 改 成 系数 精 


度 不 同 ， 使 得 对 流 层 天 顶 延 迟 估计 精度 在 不 同 高 程 范围 内 存在 一 定 的 差异 性 。 


图 5 2016-2018 年 平均 偏差 与 均 方 根 误差 在 各 高 程 范围 的 变化 
Fig.5 Variations of the bias and RMS between 2016-2018 with altitude 
3.3.2 平均 偏差 与 均 方 根 误差 的 纬度 分 布 特征 
亚洲 区 域 地 貌 和 地 形 复杂 ， 海 洋 与 陆地 分 布 比 较 广 泛 ， 气 候 复杂 多 样 ， 因 此 ZTD 变化 
也 较为 复杂 。 为 分 析 3 种 模型 的 平均 偏差 和 均 方 根 误差 在 亚洲 区 域 的 分 布 情况 , 图 6 给 
平均 偏差 和 RMS 在 各 IGS 站 的 分 布 情况 . 图 6 显示 在 亚洲 南部 和 东部 3 种 模型 的 bias 和 均 
方 根 误差 相对 较 大 ， 而 在 亚洲 南部 和 西部 误差 相对 较 小 ， 特 别 是 在 亚洲 高 纬度 地 区 
Sas+Ask+GPT3 模型 的 平均 偏差 和 均 方 根 误 差分 别 优 于 0.5cm 和 3cm， 因 此 在 亚洲 高 纬度 区 
域 和 西部 地 区 该 模型 可 以 应 用 于 高 精度 对 流 层 天 顶 延迟 研究 与 应 用 ; 在 沿海 地 区 的 IGS 站 3 
种 模型 的 误差 通 
水 汽 估算 的 精度 不 
为 进一步 分 析 3 种 模型 的 平均 偏差 和 均 方 根 误差 随 纬度 的 变化 情况 ,, 将 纬度 按 16°N ~ 


常 较 大 ， 可 能 是 沿海 地 区 易 受 海洋 气候 影响 ,水汽 丰富 且 变化 复杂 ， 导 致 对 


20°N. 20°N ~30°N, 30°N ~ 40°N, 40°N ~ 50?N 和 50°N ~ 63°N 进行 划分 ， 对 2016-2018 年 


各 纬度 范围 站 点 的 误差 进行 统计 ， 如 图 7。 图 7 显示 3 种 模型 的 均 方 根 误差 随 着 纬度 的 增加 
表现 出 减 小 的 趋势 ，Sas+Ask+GPT3 模型 在 亚洲 各 纬度 地 区 相对 于 另 两 种 模型 具有 不 同 程度 
的 改进 ， 特 别 是 在 亚洲 低 纬 度 地 区 修正 效果 更 为 显著 ， 相 对 于 SastAsk+UNB3m 和 


Sas+GPT3m 模型 ， 在 低 于 30°N 的 低 纬 度 地 区 SastAsk+GPT3 模型 的 平均 偏差 分 别 减 少 约 
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53% 和 66%， 均 方 根 误 差分 别 降 低 约 32% 和 22%; 在 高 于 30°N 的 地 区 ，Sas+Ask+GPT3 模 
型 的 平均 偏差 相对 于 另 两 种 模型 降低 约 99%， 均 方 根 误差 分 别 降 低 约 19% 和 8%。 由 图 8(a) 
还 可 以 看 到 ,在 亚洲 局 部 地 区 UNB3m 模型 的 平均 偏差 大 于 Leandro 等 在 研究 中 得 出 的 0.5cm 
的 平均 偏差 品 , 这 可 能 是 因为 UNB3m 模型 依托 北美 地 区 的 探 空 数 据 建 立 而 具有 一 定 区 域 性 
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图 6 2016-2018 年 平均 偏差 与 均 方 根 误差 在 经 纬度 上 的 分 布 


Fig.6 Distribution of the bias and RMS between 2016-2018 in longitude and latitude 


ve 一 一 一 = — ee mt 


i 


wees 
[ LAS 


图 7 2016-2018 年 平均 偏差 与 均 方 根 误差 在 各 纬度 范围 的 变化 


HE 


Fig.7 Variations of the bias and RMS between 2016-2018 with latitude 


3.4 精度 复 验 

为 进一步 验证 Sas+Ask+GPT3 模型 的 效果 及 其 时 空 变 化 特征 的 一 般 性 ,选取 位 于 中 国境 
内 的 BJFS, CHAN, CKSV, HKSL, LHAZ 和 WUH2 IGS 站 的 2019 年 1 月 1 日 至 2019 年 
12 月 31 日 的 对 流 层 天 顶 延 迟 数据 对 模型 精度 进行 分 析 评 估 , 这 6 个 IGS 站 基本 窗 盖 了 中 
的 北部 、 西 部 、 南 部 、 东 北部 、 东 南部 和 中 部 区 域 ， 具 有 典型 的 区 域 代 表 性 ， 选 取 IGS 站 
的 信息 见 表 2( 按 高 程 排序 ) ，2019 年 全 年 的 平均 偏差 和 均 方 根 误差 统计 结果 如 表 2。 

表 2 显示 CKSV 和 HKSL 沿海 区 域 的 对 流 层 天 顶 延迟 计 算 误差 较 大 ， 这 是 因为 沿海 区 
域 易 受 海洋 气候 影响 ， 空 气 潮湿 ， 对流层 湿 延 迟 的 影响 变 大 ， 经 验 模 型 对 水 汽 预 估 的 精度 不 
足 使 得 对 对 流 层 湿 延 迟 估计 出 现 较 大 误差 ，Sas+Ask+GPT3 模型 误差 明显 小 于 Sas+GPT3 模 
A, 说明 Ask 模型 在 水 汽 丰富 地 区 对 对 流 层 湿 延 迟 处 理 更 为 准确 ; Sas+Ask+UNB3m 模型 误 
差 明 显 大 于 Sas+Ask+GPT3 模型 ， 表 明 GPT3 模型 在 局 部 地 区 估算 的 气象 参数 更 为 准确 。 由 
表 2 还 可 以 看 到 , 除 在 HKWS 沿海 区 域外 , Sas+Ask+GPT3 和 Sas+GPT3 模型 的 均 方 根 误差 
随 着 高 程 的 增加 而 减 小 ，BJFS 和 CHAN 中 纬度 地 区 的 误差 小 于 低 纬度 地 区 ; 相对 于 
Sas+Ask+UNB3m 模型 ，Sas+Ask+GPT3 模型 的 平均 平均 偏差 和 均 方 根 误差 分 别 降低 了 约 
82% 和 29%， 相 对 于 Sas+GPT3 模型 ， 其 平均 偏差 和 均 方 根 误 差分 别 降 低 了 约 93% 和 14%, 
可 以 满足 中 国 区 域 对 对 流 层 天 顶 延 迟 实时 修正 的 应 用 需求 。 上 述 结果 表明 SastAsk+GPT3 
模型 的 修正 效果 显著 ， 与 前 述 结 论 相 符 。 

&2 PAE IGS 站 2019 年 年 均 平均 偏差 和 均 方 根 误差 统计 值 (单位 : cm) 


Tab.2 Statistics of yearly averaged bias and RMS in 2019 over the Chinese IGS stations (unit: cm) 


Station Latitude Longitude — SascAsk-GPT3 Sas+Ask+UNB3m Sas+GPT3 
Altitude /m 
name PN PE bias RMS bias RMS bias RMS 
WUH2 25.80 30.53 114.36 0.34 6.03 1.03 6.63 -3.26 6.81 
CKSV 59.60 23.00 120.22 -0.67 4.71 -4.79 7.61 -2.13 5.72 
HKWS 63.80 22.43 114.34 -1.01 5.00 -3.17 7.44 -4.74 7.10 
BJFS 87.41 39.61 115.89 0.19 4.39 2.81 5.68 -1.01 4.50 
CHAN 268.30 43.79 125.44 -0.28 3.70 1.37 4.88 -1.12 3.87 
LHAZ 3622.00 29.66 91.10 0.62 2.00 -1.77 4.14 1.00 2.17 
平均 -0.14 4.31 -0.75 6.06 -1.88 — 5.03 


采用 传统 对 流 层 天 顶 延 迟 改正 模型 的 建 模 思路 ， 将 GPT3 模型 和 Saastamoinen, Asknet 

模型 相 结合 ， 其 中 ，GPT3 模型 为 后 两 模型 提供 计算 对 流 层 干 延迟 和 对 流 层 湿 延迟 所 需 的 气 

RAR, 由 此 构成 Sas+Ask+GPT3 融合 模型 .选取 亚洲 地 区 66 个 IGS 站 2016-2018 年 的 GGOS 

Atmosphere fe IGS 对 流 层 天 顶 延迟 数据 对 融合 模型 进行 分 析 评 估 ， 结 果 表 明 ， 融 合 模型 在 
亚洲 区 域 对 流 层 天 顶 延 迟 修 正 效 果 显著 ， 精 度 优 于 同等 条 件 下 的 SastAsk+UNB3m 和 

Sas+GPT3 模型 ， 其 误差 具有 夏季 大 和 冬季 小 的 季节 性 特征 ， 同 时 模型 误差 随 着 IGS 站 高 程 

或 纬度 的 增加 而 减 小 ; 在 水 汽 变化 复杂 的 地 区 SastAsk+GPT3 模型 精度 明显 优 于 另 两 种 模 

型 ， 且 在 高 海拔 地 区 该 模型 的 精度 高 于 Sas+Ask+UNB3m， 说 明 GPT3 模型 比 UNB3m 模型 

在 高 海拔 地 区 估算 的 气象 参数 更 为 准确 .利用 中 国 区 域 6 个 IGS 站 2019 年 的 IGS 对 流 层 天 

顶 延 迟 数 据 对 模型 进行 进一步 验证 ， 结 果 显 示 Sas+Ask+GPT3 模型 的 误差 特征 具有 一 般 性 

本 文 研 究 表明 ， 融 合 模型 更 符合 亚洲 区 域 的 实际 情况 ,可 以 用 于 亚洲 地 区 对 流 层 天 顶 延 迟 的 

实时 修正 。 
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Abstract: This paper proposes a method fusing GPT3, Saastamoinen and askne models to 
improve the estimation accuracy of zenith tropospheric delay (ZTD) based on the conventional 
modeling strategy. In this method, Saastamoinen and askne models are used to estimate zenith 
hydrostatic and wet delay respectively, and GPT3 model is employed to provide the 
meteorological parameters including temperature, pressure, water vapor pressure, weighted 


temperature and vertical gradient of water vapor required by the two models. The ZTD data in the 


the International GNSS Service (IGS) Asia stations from 2016 to 2018 published by the Global 
Geodetic Observing System (GGOS) atmosphere and IGS are used to evaluate this proposed 
method. The results indicate that the accuracy (RMS: 4.53cm) of the Sas+Ask+GPT3 model is 
about 29% and 19% higher than that of SastAsk+UNB3m and Sas+GPT3 models respectively 
when the ZTD data from the GGOS atmosphere are served as a reference. When the IGS ZTD 
products are used as a reference, the accuracy (RMS: 4.35cm) of the model is about 25% and 14% 
higher than that of the other two models, respectively. The bias and RMS of the ZTD estimates is 
more larger in summer than in winter, showing that the error has a seasonal characteristics. 
However, the error in summer is significantly lower than the other two models. In addition, its is 
also shown that the bias and RMS decrease with the increase of altitude or latitude. Especially, the 
error in high-altitude stations is substantially lower than that of the Sas+Ask+UNB3m model. It is 
concluded that the proposed method is more potential to accurately estimate ZTD in Asia areas 
and therefore can be used for real-time ZTD correction over Asia. 

Keyword: zenith tropospheric delay; GPT3 model; UNB3m model; Saastamoinen model; Askne 


model 


